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ABSTRACT 

We present an analysis of a Chandra ACIS-S observation of the elliptical galaxy NGC 720 to verify 
the existence of a dark matter halo and to measure its ellipticity. The ACIS-S3 image reveals over 60 
point sources distributed throughout the field, most of which were undetected and therefore unaccounted 
for in previous X-ray studies. For semi-major axes a < 150" (18.2/iyg 1 kpc) the ellipticity of the diffuse 
X-ray emission is consistent with a constant value, e x « 0.15, which is systematically less than the 
values 0.2-0.3 obtained from previous ROSAT PSPC and HRI observations because of the unresolved 
point sources contaminating the ROSAT values. The Chandra data confirm the magnitude of the ~ 20° 
position angle (PA) twist discovered by ROSAT over this region. However, the twist in the Chandra 
data is more gradual and occurs at smaller a also because of the point sources contaminating the ROSAT 
values. For a > 150" out to a = 185" (22 Ah^ kpc), which is near the edge of the S3 CCD, e x and PA 
diverge from their values at smaller a. Possible origins of this behavior at the largest a are discussed. 

Overall the ellipticities and PA twist for a < 150" can be explained by the triaxial mass model of 
NGC 720 published by Romanowsky & Kochanek (which could not produce the abrupt PA twist in the 
ROSAT HRI data). Since the optical image displays no substantial isophote twisting, the X-ray PA 
twist requires a massive dark matter halo if the hot gas is in hydrostatic equilibrium. Furthermore, 
the values of e x obtained by Chandra are too large to be explained if the gravitating mass follows the 
optical light (M oc L+) irrespective of the PA twist: The M oc L* hypothesis is inconsistent with the 
Chandra ellipticities at the 96% confidence level assuming oblate symmetry and at the 98% confidence 
level for prolate symmetry. Thus, both the PA twist and the ellipticities of the Chandra image imply the 
existence of dark matter independent of the temperature profile of the gas. This geometric evidence for 
dark matter cannot be explained by alternative gravity theories such as the Modification of Newtonian 
Dynamics (MOND). 

To constrain the ellipticity of the dark matter halo we considered both oblate and prolate spheroidal 
mass models to bracket the full range of (projected) ellipticities of a triaxial ellipsoid. The dark matter 
density model, p oc (a 2 s + a 2 ) -1 , provides the best fit to the data and gives ellipticities and la errors 
of e = 0.37 ± 0.03 for oblate and e = 0.36 ± 0.02 for prolate models. Navarro-Frenk- White (NFW) and 
Hcrnquist models give similar ellipticities for the dark matter. These moderate ellipticities for the dark 
halo are inconsistent with both the nearly spherical halos predicted if the dark matter is self-interacting 
and with the highly flattened halos predicted if the dark matter is cold molecular gas. These ellipticities 
may also be too large to be explained by warm dark matter, but are consistent with galaxy-sized halos 
formed in the currently popular ACDM paradigm. 

Subject headings: X-rays: galaxies - galaxies: halos - galaxies: formation - galaxies: elliptical and 
lenticular, cD - galaxies: individual: NGC 720 - dark matter 



1. INTRODUCTION 

It is now almost 20 years that the "Cold Dark Matter" 
paradigm (CDM) , in which most of the matter in the uni- 
verse is assumed to be collisionless dark matter, has been 
regarded as the standard cosmological model of the forma- 
tion of structure in the universe. The CDM model (and 
particularly its variant ACDM) has achieved many suc- 
cesses in describing current observations, especially with 
regard to the cosmic microwave background radiation and 
the large-scale clustering properties of galaxies (e.g., Rees 
2000). Recently, however, many have drawn attention to 
the problems ACDM has describing structure on galaxy 
scales (e.g., Sellwood & Kosowsky 2002b). 



One problem with CDM is that it predicts that all galax- 
ies, including elliptical galaxies, should possess massive 
dark matter halos. But definitive evidence for such ha- 
los in isolated elliptical galaxies has proven difficult to es- 
tablish. Stellar dynamical analyses have made significant 
progress limiting the uncertainties on the radial mass pro- 
file in Es due to velocity dispersion anisotropy (e.g., de 
Zeeuw 1997; Gerhard 2002), but a recent study by Baes 
& Dejonghe (2001) argues that interstellar dust in Es ren- 
ders such mass determinations unreliable. Gravitational 
lensing studies have the potential to provide important 
constraints on the mass distributions in elliptical galaxies 
from statistical analysis of the weak shear fields of many 
galaxies (e.g., Natarajan & Refregier 2000). Statistical av- 
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eraging of the properties of strong lenses has not yet pro- 
vided clear evidence for dark matter in isolated Es (e.g., 
Keeton et al. 1998). 

There is evidence for dark matter in the outer halos 
(> 100 kpc) of elliptical galaxies obtained from kinematic 
tracers such as globular clusters, though for the few objects 
studied the inferred masses imply a surrounding group or 
cluster (e.g., in NGC 4472, Zepf 2002). Accurate con- 
straints on the radial temperature profiles of hot gas ob- 
tained by ROSAT X-ray observations have indicated sub- 
stantial amounts of dark matter for elliptical galaxies in 
the centers of groups (e.g., NGC 1399, Jones et al. 1997) 
and clusters (e.g., M87, Nulsen & Bohringer 1995). But 
because they are fainter in X-rays, ROSAT (and Einstein, 
Fabbiano 1989) did not provide temperature profiles of suf- 
ficient accuracy to clearly demonstrate the need for dark 
matter in isolated Es. 

It is clear that the evidence for dark matter in isolated 
elliptical galaxies just mentioned is not extensive, and even 
in disk galaxies the case for dark matter is not considered 
to be as strong as it once was (e.g., Evans 2002). Even for 
those galaxies where dark matter is suggested on the ba- 
sis of standard Newtonian gravity, Sellwood & Kosowsky 
(2002a) argue that all such evidence for dark matter can- 
not be distinguished from general modified gravity the- 
ories, such as the Modification of Newtonian Dynamics 
(MOND), proposed by Milgrom (1983c,a,b). 

A powerful test for dark matter that can distinguish 
dark matter from MOND in isolated Es is the "Geometric 
Test for Dark Matter" we introduced previously (Buote & 
Canizares 1994, 1996b, 1998a). The Geometric Test is a 
comparison between the ellipticitics and orientations of the 
observed X-ray isophotes with those expected if the grav- 
itating mass traces the same shape as the stellar light. 
Applying this Geometric Test to the ellipticities of the 
ROSAT data of two isolated Es and one isolated SO galaxy, 
we found strong evidence for dark matter for the E4 galaxy 
NGC 720 (99% conf.) and marginal evidence for the SO 
galaxy NGC 1332 (90% conf.) and the E4 galaxy NGC 
3923 (80%-85% conf.) - see Buote & Canizares (1998a) 
for a review. The ROSAT data of NGC 720 also displayed 
an X-ray position angle twist. But this twist could not 
be explained with triaxial dark matter models (Buote & 
Canizares 1996a; Romanowsky & Kochanek 1998) which 
cast some doubt on whether the X-ray emission indeed 
traces the potential in NGC 720. And since the evidence 
for dark matter in NGC 1332 and NGC 3923 is marginal, 
higher quality observations of these systems are required 
to assess more definitively the existence of dark matter 
halos. 

For the galaxies mentioned above we inferred from the 
ROSAT observations that the dark matter ellipticities 
range from approximately 0.4 — 0.6 with some indication 
(particularly in NGC 720) that the mass is flatter than the 
stellar distribution. These moderate ellipticities provide 
important constraints on competing models for the dark 
matter, some of which predict spherical halos (Spergel & 
Steinhardt 2000), highly flattened halos (Pfenniger et al. 
1994), and moderately flattened halos (e.g., Bullock 2002). 

4 http:/ /cxc. harvard.edu/udocs/docs/docs. html 

5 http://cxc.harvard.edu/ciao/ 

6 http://cxc.harvard.edu/cal 



Owing to the small number of measurements of the flat- 
tening of dark halos (e.g., Merrifield 2002; Sackett 1996), 
and that such previous measurements are of limited accu- 
racy, the different predictions for halo shapes have not yet 
been given a definitive test. 

We have obtained a high-resolution Chandra observa- 
tion of NGC 720 to reexamine the existence of a dark 
matter halo via the Geometric Test and to measure its 
flattening. A major issue to be addressed is the effect of 
discrete sources, such as have been found in large quan- 
tities for early-type galaxies with Chandra (e.g., Sarazin 
et al. 2000), on the conclusions drawn from lower resolu- 
tion ROSAT data. Although calculating triaxial models is 
beyond the scope of this paper, we will also measure the 
position angles of the X-ray isophotes and compare them 
to the previously published triaxial model of Romanowsky 
& Kochanek (1998). We assume a distance of 25 Mpc to 
NGC 720 for H = 70h 7Q km s" 1 Mpc" 1 (and ft = 0.3) 
in which case 1" = 0.12 kpc. 

The paper is organized as follows. We discuss the data 
preparation in §2. The treatment of point sources and the 
measurement of the X-ray ellipticitics, position angles, and 
radial profile are discussed in §3. The Geometric Test for 
dark matter is performed in §4. We examine the Chandra 
constraints on the temperature profile in §5 and the ellip- 
ticity of the dark matter in §6. Finally, in §7 we present 
our conclusions. 

2. OBSERVATIONS AND DATA REDUCTION 

NGC 720 was observed with the ACIS-S3 camera for 
« 40 ks during AO-1 as part of the guaranteed-time- 
observer program (PI Garmire). The events list was cor- 
rected for charge-transfer inefficiency according to Towns- 
ley et al. (2002), and only events characterized by the stan- 
dard ASCA grades 4 were used. The standard CIAO 5 soft- 
ware (version 2.2.1) was used for most of the subsequent 
data preparation. 

Since the diffuse X-ray emission of NGC 720 fills the en- 
tire S3 chip, we used the standard background templates 6 
to model the background. When attempting to use the 
LC_CLEAN script to clean the source events list of flares 
with the same screening criteria as the background tem- 
plates, only 12 ks of data remained. Examination of the 
original light curve reveals only one significant flare (~ 6 ks 
worth of data) near the beginning of the observation, but 
the light curve does gently rise and fall during the rest of 
the observation so that most of the data does not lie within 
the standard 3-sigma clip. Consequently, we excluded only 
the obvious flare leaving a total exposure of 34.4 ks. 

Because the quiescent background can vary typically by 
~ 10% between observations, and we do not use exactly 
the same screening criteria as the templates, it is necessary 
to normalize the background templates to match the NGC 
720 observation. Since NGC 720 has a gas temperature 
~ 0.6 keV there is little emission from hot gas for energies 
> 5 keV. Consequently, we normalized the S3 background 
template by comparing source and background counts in 
the 5-10 keV band extracted from regions near the edges 
of the S3 field. 
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We binned the cleaned S3 events lists for the source 
and background into images of 1" x 1" pixels (actually 
0.984" x 0.984" pixels, but we shall round off throughout 
the paper). Since background dominates at high energies 
we selected events with energies 0.3-3 keV. Likely prob- 
lems with the spectral calibration near the low end of the 
bandpass (see §5) are not important for the imaging anal- 
ysis. 

Although vignetting is a small effect within « 3' of the 
aim point of the S3 where our analysis is concentrated, we 
explored the viability of using an exposure map to flat- 
ten the image. The inherent difficulty in constructing a 
single exposure map for the image of a diffuse source for 
instruments like AC1S has been discussed by Davis (2001) 
and Houck 7 . Since, however, the diffuse gas of NGC 720 is 
remarkably isothermal (see §5) we are able to construct a 
well-defined "weighted" exposure map following the pro- 
cedure on the CIAO website. (Note we find no tangible 
differences in any results mentioned below if instead we 
use a monochromatic exposure map evaluated at 0.77 keV 
which is the peak of the counts spectrum.) 

Unfortunately, the current method for flat-fielding im- 
ages introduces spurious features along node boundaries 
where the intrinsic exposure is small relative to the rest 
of the CCD. When the exposure map is divided into the 
raw image to create an exposure-corrected image, regions 
of relatively low exposure (e.g., node boundaries, bad 
columns, and edges) require large corrections. Any in- 
trinsic errors in the exposure map in these regions, and 
any noise in the data, will be amplified as a result. 

For NGC 720 we find that dividing by the exposure map 
amplifies the linear distortion along the node 1-2 bound- 
ary that reaches within w 30" of the center of NGC 720. 
Ellipticities and position angles computed between radii of 
- 30" -60" are systematically lower by 20% -30% with re- 
spect to their values obtained from the un-flattened image. 
We emphasize that in the raw image the linear distortion 
is very faint and very difficult to see, and imperceptable in 
the smoothed image displayed in Figure 1; i.e., this feature 
has negliblc effect on the ellipticities and position angles 
considering the statistical uncertainties on these quanti- 
ties (§3.2.2). Consequently, since we have a case where 
"the cure is worse than the disease" , when computing el- 
lipticities and position angles we do not use the exposure- 
corrected image. But since this non-circular, low-surface 
brightness feature amplified by the exposure map has no 
perceptible effect on the azimuthally averaged radial pro- 
file, we do use the exposure-corrected image to account 
for the small effects of (large-scale) telescopic vignetting 
for analysis of the image radial profile. 

3. IMAGE ANALYSIS 

In Figure 1 we display the adaptively smoothed S3 image 
of NGC 720 within a radius of 150" of the galaxy center 
with some isophotal contours overlaid. (This smoothing is 
done only for display purposes. All analysis discussed be- 
low is performed on the unsmoothed image; cf. Figure 1 of 
Jeltema et al. 2002). The most striking difference between 
the Chandra and ROSAT images (e.g., Buotc & Canizares 
1998a) is the large number of point sources clearly detected 
by Chandra. Many of these sources are located within 1' 

7 see http://cxc.harvard.edu/ciao/documents_manuals.html 



of the galaxy center. 

The morphology of the diffuse emission within ~ 100" of 
the galaxy center is similar to that described for previous 
ROSAT observations. The isophotes are moderately flat- 
tened and the orientations of the outer isophotes appear 
to be tilted with respect to the innermost regions. But 
caution must be exercised when viewing this adaptively 
smoothed image since small-scale features in the diffuse 
emission have low statistical significance (especially ap- 
parent wiggles in the isophotes). Below we shall quantify 
the morphology of the diffuse emission using moments of 
the image. 

3.1. Removal of Point Sources 

Before we can analyze the properties of the diffuse emis- 
sion we must remove the point sources. (Analysis of the 
properties of these sources is presented in a companion pa- 
per - Jeltema et al. 2002.) However, it is insufficient only 
to remove the sources since holes in the image still affect 
significantly quantities like the ellipticity that are derived 
from second (or higher) moments of the image. To obtain 
reliable measurements of the ellipticity and position angle 
throughout the source-free image it is necessary to replace 
each embedded source with a faithful representation of the 
local diffuse X-ray emission surrounding the source. 

We located point sources with the CIAO task WAVDE- 
tect (a wavelet detection algorithm) using the default 
parameter settings. A total of 64 sources were detected 
over the entire S3 field excluding the center of NGC 720. 
These include the 41 sources within 2.5' of the galaxy cen- 
ter discussed by Jeltema et al. (2002). 

To remove these sources and "fill" in the holes with dif- 
fuse emission we used the CIAO task DMFILTH. This task 
takes as input source and background region definitions. 
The counts within the source regions are replaced with val- 
ues sampled from the background regions. Within ~ 3' of 
the galaxy center we defined the source extraction regions 
as circles of radius 2". The corresponding concentric back- 
ground regions all have radius 4". The average background 
counts in each region was determined by the polynomial 
method in DMFILTH. Since in several cases the background 
region of one source overlapped another (especially within 
~ 30" of the galaxy center) we iterated runs of dmfilth 
using the cleaned images of previous runs as input to sub- 
sequent runs. We found that the results converged after 
only three runs. 

For radii larger than ~ 30" we found that dmfilth 
performed well. In this context "well" means that the 
measured ellipticity and position angle are smoothly vary- 
ing functions of radius; i.e., at the radius of a filled- in 
source the ellipticity and position angle do not deviate 
substantially from values obtained at adjacent source-free 
radii. Unfortunately, for radii < 30" we found that the 
ellipticity and position angle did jump "discontinuously" 
near some filled-in sources, especially at the smallest radii 
< 15" where source background-region overlap (as noted 
above) is most pronounced. (Note that such jumps, by 
accident, are not apparent for the ellipticity and position 
angle at the particular radii < 30" plotted below.) 

At small radii the relative effect of the sources on de- 
rived surface brightness parameters is necessarily greater 



Fig. 1. — (Left) False-color ACIS-S3 image of NGC 720 in the 0.3-3 keV band. The image has been adaptively smoothed using the CIAO task 
CSMOOTH with default parameters. No exposure-map correction or background subtraction has been applied. Contours are spaced according 
to the square root of the intensity. (Right) Same image except the point sources have been removed. The contour levels are the same as the 
original image. Each image is oriented so that celestial N is up and E is to the left. 



because the area of the source becomes increasingly more 
significant relative to r 2 . It is likely that at such small 
radii the approach of simply filling in the source region 
with a spatially uniform count distribution does not faith- 
fully (enough) represent the shape the local diffuse emis- 
sion. We defer consideration of more sophisticated source- 
replacement algorithms to a future study. In our analyses 
below we treat the region r < 30" with caution. 

As should be expected, residual source-replacement ef- 
fects in these inner regions are much less important for the 
radial profile. We shall therefore use the whole image for 
construction and analysis of the radial profile. 

3.2. e x and PA 

3.2.1. Method 

The flattening of the X-ray isophotes is of central im- 
portance to our study of the flattening and concentration 
of the total gravitating matter in NGC 720. The standard 
parameter used to denote the flattening of an individual 
isophote is the ellipticity, 1 — b/a, where a is the semi-major 
axis and b is the semi-minor axis. However, fitting perfect 
ellipses to the X-ray isophotes is not necessarily justified 
since the isopotential surfaces generated by an elliptical 
mass distribution are not perfect ellipsoids. In addition, 
over most of the X-ray image the surface brightness of the 
diffuse gas is < 1 counts per pixel. This condition requires 
that relatively large areas (in excess of a single-pixel width 
isophote) must be averaged over to obtain interesting con- 
straints on the image flattening. 

We quantify the image flattening using the method de- 
scribed by Carter & Metcalfe (1980) and implemented in 
our previous study of the ROSAT image of NGC 720 (e.g. 
Buote & Canizares 1994). This iterative method is equiv- 
alent to computing the (two-dimensional) principal mo- 
ments of inertia within an elliptical region. The ellipticity, 
e x , is defined by the square root of the ratio of the princi- 
pal moments, and the position angle, PA, is defined by the 
orientation of the larger principal moment. If these mo- 
ments are computed within an elliptical region where the 
image is perfectly elliptical with constant ellipticity and 



orientation, then e x represents a true ellipticity and PA is 
the true orientation of the major axis. If the image is not 
perfectly elliptical within this region, then e x and PA are 
average values weighted heavily toward the edge of the re- 
gion; i.e., e x provides a useful measure of image flattening 
which does not assume the image to be perfectly elliptical. 

Following our previous studies we compute e x (a) using 
all image pixels interior to the ellipse defined at a (obtained 
via iteration). We did investigate using elliptical annuli to 
provide a more direct measurement of the variation of e x 
with a. Unfortunately, because of the low counts per pixel 
noted above the radial fluctuations of e x , PA, and centroid 
were considerably larger when using elliptical annuli. 

We estimate uncertainties on e x and PA using a Monte 
Carlo procedure. The counts in each pixel are randomized 
assuming poisson statistics. Then e x and PA are computed 
precisely as done for the original image. After performing 
100 such realizations and corresponding "measurements" 
of the image we compute the standard deviations of e x and 
PA at each a for the 100 runs. We take these standard de- 
viations to be the la errors. 

We have computed e x and PA out to a = 185" from 
the Chandra image such that the background-subtracted 
counts in each aperture increase by « 500 for each a. (The 
edge of the S3 chip corresponds to r w 215".) Because we 
use all image pixels interior to a to compute e x and PA, 
their values computed for a given a are correlated with 
values computed at adjacent inner a. This means that the 
error bars on e x and PA for adjacent a are not fully inde- 
pendent. But we emphasize that these quantities, being 
derived from second moments of the image, are weighted 
heavily by the pixels near a. 

3.2.2. Results 

In Figure 2 we show e x and PA computed from the Chan- 
dra image before and after removing the point sources. 
(Note the error bars for values computed from the image 
with sources are not shown for clarity but are of similar 
magnitude to those obtained from the source-free image.) 
For the image with point sources, within a < 80" e x varies 
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Fig. 2. — (Left panel) e x as a function of semi-major axis computed from the Chandra image with point source removed (circles and error 
bars; blue). The (red) boxes are the values of e x obtained from the image which includes the point sources. The (green) dotted error bars 
are the values obtained with the ROSAT HRI (Buote & Canizares 1996a) (Right panel) PA as a function of a computed from the Chandra 
image with point sources removed (circles and error bars; blue). The (red) boxes are the values of e x obtained from the image which includes 
the point sources. Note that because we use all image pixels interior to a to compute e x and PA, the values and the error bars for adjacent a 
are not independent. We express a in kpc on the top axis. 
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Fig. 3. — (Left paneJ) e x as a function of semi-major axis computed from the source-free Chandra image (circles and error bars; blue). The 
optical B-band cllipticity is indicated by the (green) dashed line (see §4.2). The values of e x expected for a model where the gravitating mass 
follows the optical light is shown by the (red) squares; i.e., the oblate M oc L* model mentioned in the text. (Right panel) PA as a function 
of a computed from the source-free Chandra image (circles and error bars; blue). The optical B-band value (green) is also shown. We express 
a in kpc on the top axis. 



wildy with a then settles down to e x w 0.2 until it declines 
for a > 150". This erratic behavior is a stark constrast to 
the slowly varying values of e x computed from the source- 
free image. We interpret this result as a strong affirmation 
of our success at removing and replacing the point sources 
over this range in a for the purpose of computing e x . 

The PA generally does not exhibit such dramatic varia- 
tions with a in either case, though the PA values computed 
from the image with sources (for the most part) systemati- 
cally exceed those computed from the source- free image by 
~ 10°. This relative similarity between the PAs computed 
from the raw and source-free images coupled with the lack 
of strong fluctuations with a suggest that PA is less sensi- 
tive than is e x to contaminating point sources. However, 
with the present data we cannot rule out the possibility 
that PA is more sensitive than e x to contaminating point 



sources and is therefore only slightly affected by our at- 
tempt to exclude and replace the sources. Since the values 
of e x and PA are necessarily interwined, and we are confi- 
dent of the success in computing e x , we do not believe the 
PA values computed from the source-free image are still 
substantially affected by embedded point sources. Never- 
theless, to be conservative we shall address below the need 
for dark matter in NGC 720 using e x and PA separately. 

Let us focus on the source- free image (Figure 3). As ex- 
plained above in §3.1 the best-fitting values of e x and PA 
computed within a < 30" must be considered tentative be- 
cause of residual source-replacement errors. Nevertheless, 
certain trends in these data appear to be robust. First, 
for small a we have, e x « 0.2 — 0.3, consistent with the 
optical isophotes in this region within the relatively large 
estimated errors. Second, the position angle is consistent 
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with the optical value, PA « 142° (N-E), for a < 30". 

For 30" < a < 150" the ellipticity is consistent with a 
constant value, e x w 0.15. This value is notably smaller 
than the best-fitting values of « 0.20 — 0.30 obtained with 
the ROSAT PSPC and HRI (Buote & Canizares 1994, 
1996a, 1998a) between 30" < a < 110"; the HRI values 
are reproduced in Figure 2. The overall consistency of the 
e x values obtained from the ROSAT image and Chandra 
image with sources demonstrates that the embedded point 
sources are responsible for the ROSAT values. The most 
important are the group of point sources located « 70" to 
the NW of the galaxy center. These sources were not iden- 
tified and removed from the ROSAT images (and evaded 
detection from our symmetry tests because they are lo- 
cated near the X-ray major axis). Note that the source- 
free Chandra e x values are smaller than obtained from the 
HRI for 50" < a < 70". This discrepancy is not resolved 
by simply degrading the Chandra data by smoothing it 
with the HRI PSF. We have not definitively identified the 
origin of this discrepancy, but we suspect it is assocated 
with the different flux sensitivities of the HRI and ACIS; 
i.e., flat spectrum point sources will be more prominent 
with respect to the diffuse hot gas in the ACIS data. 

The magnitude of the PA twist discovered by ROSAT is 
confirmed by the source-free Chandra data, but the pre- 
cise radial variation of PA is different because of the un- 
resolved point sources present in the ROSAT data. For 
80" < a < 150" the position angles, PA w 110°, measured 
by Chandra are consistent with the ROSAT values. How- 
ever, for a - 60" the ROSAT HRI PA jumps to a value 
PA « 140° whereas the source-free Chandra PA increases 
more gradually from PA 120° to PA « 135° for decreas- 
ing a within a < 60". The similarity of the PA values 
obtained for the Chandra image with sources and the HRI 
demonstrate that the embedded point sources are respon- 
sible for the inflated PA values. The key point sources for 
these inflated PA values are located ~ 45" to the SW of 
the galaxy center along the minor axis (see Figure 1 above 
and Figure 1 of Jcltema et al. 2002). 

At large radius (150" < a < 180") e x increases steadily 
to a value of w 0.20 while the position angle decreases to 
a value of w 80° . These changes are highly significant sta- 
tistically and are obvious upon inspection of the smoothed 
image displayed in Figure 1. For a > 120" the PSPC and 
Chandra values of e x agree within errors. This suggests 
that the effect of unresolved sources on e x is less impor- 
tant at large radius as might be expected since the PSF 
FWHM compared to o decreases with increasing a. How- 
ever, the PA measured by the PSPC at large radius does 
not appear to twist so much as seen in the Chandra data: 
PSPC PA values are 102° ± 9° at a = 150" and 107° ± 5° 
at a = 225" (Buote & Canizares 1994). 

The origin of this Chandra-ROSAT PA discrepancy at 
large radius is unclear, but we would tend to believe the 
Chandra values because of the issue of point sources in 
the ROSAT data. It must be acknowledged, however, that 
since the edge of the S3 chip is at a « 215" there could be 
problems with flat-fielding which we think unlikely based 
on our own investigations with exposure maps in this re- 
gion. It is probably necessary to use the wide-field XMM 
EPIC CCDs to resolve this issue. 



3.3. Radial Profile 

The radial profile (£ x (r)) is the final property derived 
from the X-ray image that is required to constrain the 
shape of the gravitating matter distribution. We define 
the image radial profile by first summing the background- 
subtracted counts in concentric circular annuli with radii 
determined by the semi-major axes used to compute e x 
and PA (§3.2.2). Then dividing the counts in each annulus 
by the area of the annulus we arrive at the radial profile. 
The center of each annulus is the same and is taken to be 
the centroid computed within a circle of radius 20" placed 
initially at X-ray peak. (The centroid varies by < 4" for 
r < 185".) There are « 12000 background-subtracted 
counts within r = 185". 

We plot S x (r) in Figure 4. The radial profile is quite 
smooth over the entire range investigated. The some- 
what bumpy HRI profile (Buote & Canizares 1996a) within 
r < 30" can be attributed to noise and to the point sources 
detected in that region by Chandra. Fitting the Chandra 
profile with a (3 model yields r c ~ 4.6" and (3 « 0.42 in 
agreement with the HRI values. The [3 model is an excel- 
lent qualitative (i.e., visual) fit to S x (r), but because the 
error bars are very small the formal quality of the fit is 
rather poor: \ 2 = 44.3 for 16 degrees of freedom (dof). 
The most important fit residuals occur in the inner bins. 
These residuals are small when expressed as a ratio of data 
to model (i.e., < 10%), but are relatively large when ex- 
pressed in terms of a x 2 difference between data model. 
If, for example, the inner 3 bins are excluded (r < 16") we 
obtain x 2 = 24.9 for 13 dof. 

4. GEOMETRIC TEST FOR DARK MATTER 

4.1. Basic Idea and Assumptions 

Previously, we have shown that it is possible to test the 
hypothesis that the gravitating mass traces the stellar light 
using only the X-ray imaging data of the diffuse hot gas 
- no knowledge of the value of the gas temperature or a 
possible spatial temperature gradient is required (Buote & 
Canizares 1994, 1996b, 1998a). The key assumptions are 
that the hot gas is single-phase and in hydrostatic equilib- 
rium. If these conditions hold, then the gas volume emis- 
sivity, j x , and the potential of the gravitating matter, <£, 
must have identical shapes in three-dimensional space ( "X- 
Ray Shape Theorem"); i.e., in 3D the surfaces of constant 
j x are surfaces of constant $ independent of the temper- 
ature profile of the gas. Therefore, one may test the hy- 
pothesis that the 3D gravitating mass is distributed in the 
same way as the 3D stellar light, L*, obtained from depro- 
jecting the optical image. This "Geometric Test for Dark 
Matter" compares the shapes of the isopotential surfaces 
computed assuming M oc L* to the shapes of the isoemis- 
sivity surfaces of j x , obtained from deprojecting the X-ray 
image. 

It is important to remember that not only are the shapes 
of the isopotential (and therefore isoemissivity) surfaces 
sensitive to flattening of the gravitating mass, but they 
also reflect its radial distribution. This is illustrated by 
considering a multipole expansion of the potential for an 
arbitrary mass distribution of finite extent. Near the cen- 
ter of the mass many high-order multipole terms may con- 
tribute resulting in a highly flattened potential in that 
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region. As one moves away from the center the higher 
order terms must decay and eventually (and usually very 
rapidly) give way to the spherical monopole term. Hence, 
an ellipticity gradient in the potential, and therefore j x , 
also probes the radial distribution of gravitating mass. 

This Geometric Test for dark matter that is distributed 
differently from L± is also a powerful test of alternative 
gravity theories such as MOND. The fundamental MOND 
equation is n.g M = .9n, where <?n is the Newtonian force 
vector, (Jm is the MOND force vector, and /i is a scalar 
function depending only on the ratio gm/ao ~ the magni- 
tude of the MOND force vector over a constant. Conse- 
quently, for a given mass distribution the shapes of isopo- 
tential surfaces in MOND are the same as the shapes of 
isopotential surfaces in Newtonian theory. So if the Geo- 
metric Test implies dark matter using Newtonian theory it 
follows that dark matter is also required using the MOND 
theory. See Buote & Canizares (1994) for a more extensive 
discussion of the Geometric Test, MOND, and caveats. 

The assumption of a single-phase gas for NGC 720 ap- 
pears to be well justified. First, the Chandra spectrum (see 
§5) is consistent with an isothermal plasma. Second, re- 
cent observations with Chandra and XMM have found that 
even the gas in cooling flow groups and clusters, systems 
that might have been expected to be highly multiphase, 
is consistent with a single-phase medium (e.g., Bohringer 
ct al. 2001; Buote ct al. 2002). 

There are also several lines of evidence which suggest 
that the approximation of hydrostatic equilibrium is a 
good one (see also §7.4). Although it is natural to won- 
der whether the gas is disturbed by tidal or ram pressure 
effects, it must be remembered that NGC 720 is quite iso- 
lated from other large galaxies (Dressier et al. 1986). Aside 
from the PA twist there are no indications of external in- 
fluences in the diffuse X-ray emission; e.g., centroids com- 
puted for the elliptical apertures in §3.2.2 are all consistent 
within 3" - 4". 

We should also emphasize that the wiggles and other 
irregular features displayed by the smoothed X-ray 
isophotes in Figure 1 are not statistically significant and 
thus do not necessarily imply departures of the gas from 
hydrostatic equilibrium. However, even if those wiggles 
do in fact arise from real non-equilibrium motions in the 
gas, the hydrostatic assumption is still very likely a good 
one. Previously, in Buote & Tsai (1995) we analyzed the 
X-ray isophotes of a galaxy cluster formed in an N-body 
/ hydrodynamical simulation. It was found that relatively 
soon after major mergers (i.e., 1-2 crossing times) the clus- 
ter settled down to a quasi-relaxed state where the hot 
gas traced approximately the shape of the gravitational 
potential even though the shapes of the X-ray isophotes 
displayed prominent irregular features. 

As in our previous studies we shall assume negligible 
rotation of the hot gas. We expect this to be a good ap- 
proximation for several reasons. Most importantly, like 
most massive elliptical galaxies there is negligible stel- 
lar rotation; NGC 720 has a mass-weighted stellar rota- 
tional velocity of 35 km s" 1 (Busarello et al. 1992) and 
v/cr* = 0.15 (Fried & Illingworth 1994). And since much, 
perhaps most, of the hot gas is produced by stellar mass 
loss (Mathews 1990), the hot gas should also rotate neg- 
ligibly. In principle, if the hot gas flows inwards while 



conserving angular momentum then the resulting spin-up 
could translate into dynamically important rotational ve- 
locities. (This would also apply to any primordial gas flow- 
ing in from large radius.) Both theory and the observed 
lack of highly flattened X-ray images of galaxies and cen- 
trally E-dominated groups indicate that it is highly un- 
likely that angular momentum is conserved in the hot gas 
(Nulsen et al. 1984; Brighenti & Mathews 2000). At any 
rate, the lack of evidence for cooling flows noted above 
argues against inflowing gas. It would also be difficult to 
reconcile the X-ray PA twist in NGC 720 with substantial 
solid body rotation around a single symmetry axis. Thus, 
all evidence suggests negligible rotation of the hot gas. 

Finally, we shall neglect the self-gravity of the hot gas in 
our calculations. We have shown previously using ROSAT 
data (Buote & Canizares 1994) that the mass of the hot 
gas is < 1% of the gravitating mass over the entire region 
studied. 

4.2. Current Objective and Implementation 

The X-ray PA twist provides immediate geometric ev- 
idence for dark matter since the stellar isophotes do not 
display any substantial twist. However, following our dis- 
cussion in §3.2.2 regarding the removal and replacement 
of point sources on the derived PA values, here we investi- 
gate whether the ellipticities and radial profile of the X-ray 
image corroborate the evidence for dark matter provided 
by the PA twist. We will focus our analysis of e x in the 
region 35" < a < 139". The lower limit is chosen to avoid 
residual point source-replacement issues (§3.1). The up- 
per limit avoids the outer low-surface brightness regions 
where e x and PA diverge from their values at smaller o. 
(The origin of this divergence is unclear. Some possibili- 
ties are discussed in §7.) Since S x (r) is quite insensitive 
to residual source-replacement issues and small variations 
in e x and PA, we will use its constraints over the entire 
range r < 185". 

The first ingredient required for the Geometric Test 
is the potential computed assuming M oc L t . Us- 
ing the effective radius of R e = 52" obtained from op- 
tical surface photometry of NGC 720 (Burstein ct al. 
1987), we adopt a corresponding effective semi- major axis, 
a e = R e / \J\ — (e op t), where (e op t) = 0.31 is the intensity 
weighted radial average of the optical B-band isophotal el- 
lipticities, e pt, computed from the tables of Peletier et al. 
(1990) and Lauer et al. (1995) (see Figure 3). We approx- 
imate with a spheroidal Hernquist model (Hernquist 
1990, 1992) with scale length, a* = a e /1.8153 = 34.5", 
and constant ellipticity, = 0.31. 

We will consider both oblate and prolate spheroids to 
bracket the range of projected ellipticities of a triaxial el- 
lipsoid. With no loss in the robustness of the Geometric 
Test we shall also assume that the symmetry axis lies in 
the plane of the sky. According to the "X-Ray Shape The- 
orem" the unknown inclination angle i must be the same 
for and j x if M oc L+. Since an inclined spheroid is 
necessarily rounder in projection than when viewed edge- 
on, deprojection assuming i = 90° will yield and j x 
with rounder shapes than in actuality. This effect tends 
to smear out differences in the shapes of and j x and 
therefore makes it harder to falsify the Mal t hypothe- 
sis; i.e., the assumption of i — 90° is conservative in this 
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context. 

After calculating the next step is to compare the 
shapes of the isopotential surfaces to the shapes of the 
isoemissivity surfaces of j x . To perform this compari- 
son, we parameterize j x with a f3 model 8 along the major 
axis, z (i.e., the long axis in the sky plane): j x (0,0, z) oc 
(r c 2 + z 2 )~ 3/3 where r c is the "core radius". If r c and (3 
are given, then any j x (x, y, z) can be computed by identi- 
fying the surfaces of constant j x with those of i.e., for 
any point (x, y, z) we find the corresponding value Z such 
that &+(x,y, z) = &+(0,0, Z). Then we set j x (x,y,z) = 
j x (0,0, Z) at which point we add a 'V to j x (i.e., j*) to 
emphasize that its isoemissivity surfaces have been aligned 
with 

We obtain the appropriate values of r c and (3 by pro- 
jecting j* along the line of sight to obtain £*. Then 
we compute S*(r) for this M cx L± model analogously 
to the Chandra data and perform a x 2 fit to the data. 
Interestingly, this flattened (3 model is a slightly bet- 
ter fit (x 2 = 33.9) than the spherical version reported 
in §3.3. The best-fitting parameters and lcr errors are, 
r c = 5.3" ± 0.3" and (3 = 0.414 ± 0.002. These lcr errors 
are computed using a Monte Carlo procedure. We random- 
ize the best-fitting S*(r) according to Gaussian statistics 
and fit the model as described above. After twenty such 
realizations are performed we compute the standard de- 
viations of r c and (3 obtained from the 20 runs. We take 
these standard deviations to be the lcr errors. 

The X-ray image ellipticities for this M oc model, 
denoted by e x , are also computed from S* analogously to 
the Chandra data. We also compute the lcr errors on e* 
using the 20 Monte Carlo runs mentioned above. For the 
oblate M oc L* model we plot e* and the lcr errors in Fig- 
ure 3 indicated by the "M oc L sta rs Prediction" . For most 
values of a the lcr errors are smaller than the size of the 
value marker. 

Overall the values of e x are smaller than observed. For 
a < 50", e* usually lies within the error bars estimated for 
e x computed from the Chandra image. For larger values 
of a, e* increasingly deviates from the observed e x val- 
ues; i.e., the deviation of M oc prediction from the 
observations becomes pronounced outside of <~ R e . The 
primary reason why the M oc L± model cannot produce 
the X-ray ellipticities is because L* is too centrally concen- 
trated. Consequently, the spherically symmetric monopole 
term dominates the stellar potential and predicts values of 
e x that are smaller than observed, particularly for a > 60". 
(Note the downturn in e* at small a is an artifact of the it- 
erative moment method used to calculate the ellipticities. 
The isophotal ellipticities of the M oc L* model continue 
to increase toward the center.) 

To quantify the significance of this deviation between e* 
and e x we perform a \ 2 test defined in the usual way, 



particular a are correlated with values at adjacent inner a, 
we cannot interpret the result with the standard x 2 null 
hypothesis probability which assumes uncorrelated errors 
for each point. 

Consequently, we construct our own \ 2 probability func- 
tion that accounts for the special correlations associated 
with our analysis. To do this we compute x 2 m a world 
where S* is the real X-ray image of NGC 720. That is, 
starting with S* generated using the best fitting values 
r c and (3 we compute e*; these are just the best- fitting 
values (e*)i shown in Figure 3. Then we follow the same 
procedure used to compute the uncertainties on e x from 
the real data (§3.2.1); i.e., the model counts in each X* 
image pixel are randomized assuming poission statistics, 
and then e* is computed in the same manner as above. 
Hence, for each Monte Carlo realization of £* we obtain 
ellipticities (e*)^, where i indicates the value of a and j 
indicates the particular realization. 

We perform 10 4 such Monte Carlo realizations. Using all 
of these runs we compute the standard deviations for each 
(e*)^. (The magnitudes of these lcr errors are very simi- 
lar to those obtained for the actual data shown in Figure 
3.) Now for each realization j we calculate {x 2 )j replacing 
(e x )i and {a x )i in equation (1) with respectively (e*)^ and 
their associated standard deviations. This set of 10 4 values 
of {x 2 )j defines the special x 2 probability distribution. 



4.3. Results 

The x 2 test reveals that the oblate MocL 4 model is in- 
consistent with e x and S x (r) computed from the Chandra 
image at the 96% confidence level. The prolate M oc L± 
model is inconsistent at the 98% confidence level. These 
results correspond to e x values computed within the re- 
gion 35" < a < 139" as noted at the beginning of this 
section. Since the deviation between the M oc L± model 
and data increases with a (Figure 3), the significance of 
the discrepancy increases if we exclude data at smaller 
a; e.g., if we consider instead e x values computed within 
the region 66" < a < 139", the oblate (prolate) M oc L* 
model is inconsistent with the Chandra data at the 98% 
(99%) confidence level. (If all e x values computed within 
a = 185" are used then the magnitude of the discrepancy 
is > 99.99% for both oblate and prolate models.) 

Thus, using only the ellipticities to quantify a discrep- 
ancy, the Geometric Test demonstrates that the M oc L± 
hypothesis is unsatisfactory and that dark matter is re- 
quired - independent of the temperature profile of the hot 
gas. This corroborates the evidence for dark matter pro- 
vided by the PA twist. As discussed above, such "Geo- 
metric" evidence for dark matter cannot be explained by 
general modified gravity theories such as MOND. 

5. TEMPERATURE PROFILE 

Unlike the Geometric Test, to construct detailed mod- 
els of the gravitating matter distribution we must have 
some knowledge of the temperature profile of the gas. We 
extracted the ACIS spectra in three circular annuli con- 
taining w 3000 background-subtracted counts over 0.3-3 
keV: 0" - 30", 30" - 70", and 70" - 140". We fit these 

8 Here the (3 model is used only as a convenient and accurate parameterization of S x (r). We do not attribute the physical significance to the 
model (or make many of the assumptions) that inspired its creation by Cavaliere & Fusco-Femiano (1978). 
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where i represents each value of a included in the sum, cr* 
is the standard deviation for e* mentioned above, and cr x 
is the standard deviation associated with e x (§3.2.1). But 
since values (and standard deviations) of e* and e x for a 
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spectra with an APEC thermal plasma 9 modified by Galac- 



tic absorption (N§^ = 1.55 x 10 



20 



J ). The spectral 



& Canizares 1994, 1996b, 1998a, b) based on the pioneer- 
ing approach of Binney & Strimpel (1978). We solve the 
equation of hydrostatic equilibrium for the gas density as- 
suming a single-phase, isothermal, ideal gas, 
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fitting was performed with xspec vll.l.Ov (Arnaud 1996). 

Attempts to fit this model to the 0.3-3 keV data were 
unsuccessful for two reasons. First, below ~ 0.5 keV the 
data lie well below the model and can only be reproduced if 
absorption column densities of iV^ al w 10 x 10 20 cm~ 2 are 
fitted. This effect is undoubtedly related to existing prob- 
lems in the low-energy calibration mentioned on the Chan- 
dra web site. We suspect that because the diffuse gas of 
NGC 720 has a relatively low temperature (~ 0.6 keV) this 
calibration problem is more serious for NGC 720 than hot- 
ter systems. And it is more serious for the diffuse gas than 
for the harder discrete sources (see Jeltema et al. 2002). 

If we exclude data below 0.7 keV then much better fits 
can be obtained for Galactic absorption. However, there is 
also excess emission above the model for energies > 2 keV. 
Since this excess is most pronounced in the outer annulus 
it is very likely associated with the background normaliza- 
tion. (Recall that the background varied throughout the 
observation - §2.) To largely avoid this excess we then 
excluded data above 2 keV. 

Fitting the models over 0.7-2 keV yields temperatures 
of ~ 0.6 keV in each annulus. Since, however, we have 
restricted the bandwidth to such a narrow region (which 
also does not enclose the energy corresponding to the best- 
fitting temperature) the constraints on the spectral model 
are relatively weak. In particular, the best-fitting metal- 
licity is ~ O.IZq but it is not well constrained and depends 
sensitively on the chosen upper limit of the bandpass. The 
temperatures have best-fitting values of 0.62, 0.58, and 
0.62 keV respectively in the three annuli; each has la er- 
rors of « 0.1 keV. These values are fully consistent with 
the ROSAT and ASCA values obtained in these regions 
(Buote & Canizares 1994, 1997; Buote & Fabian 1998) 

Hence, the Chandra data confirm that the hot gas in 
NGC 720 is consistent with being isothermal. Given the 
extenuating issues mentioned above, this is the only robust 
conclusion we believe can be drawn from the spectrum of 
the diffuse gas at this time. 

6. ELLIPTICITY OF THE DARK MATTER HALO 

The Geometric Test established the need for a dark mat- 
ter halo that is distributed differently from the stars. Now 
we wish to find the range of halo ellipticities that are con- 
sistent with the Chandra data. Since we are concerned 
primarily with the flattening of the halo we shall follow 
the approach adopted in §4 and consider both oblate and 
prolate spheroids to bracket the range of projected ellip- 
ticities of a triaxial ellipsoid. We defer consideration of 
the more computationally expensive triaxial models to a 
future study, but in §7 we discuss briefly the suitability of 
a previously published triaxial model for NGC 720. 

6.1. Preliminaries 

To constrain the halo ellipticities we follow the proce- 
dure outlined in our previous related studies (e.g., Buote 

9 http://hea-www.harvard.edu/APEC/ 

10 In principle the method used to test the MocL t model in §4 can be applied to any mass model to constrain its cllipticity and scale length. 
This approach has the advantages of not requiring one to specify (1) the temperature profile and (2) how the gas density p g is disentangled 
from S x . The primary disadvantage is that the method based on the Geometric Test does not, for a general mass model, constrain the mass 
scale length as well as the approach adopted in this section. 



where T = —pm p $(0)/kBT and where, e.g., p g = 
p g (x)/p g (0). The assumption of an isothermal gas is an 
excellent approximation for obtaining constraints on the 
halo cllipticity of NGC 720 10 . It has been shown previ- 
ously by several studies (Strimpel & Binney 1979; Fabri- 
cant et al. 1984; Buote & Canizares 1994; Buote & Tsai 
1995) that constraints on the cllipticity of the gravitat- 
ing matter assuming an isothermal gas are very accurate 
even if there is a sizable temperature gradient. (This is a 
direct consequence of the X-ray Shape Theorem noted in 
§4.) In particular, for NGC 720 the Chandra spectral con- 
straints (§5) indicate a temperature profile that cannot 
deviate substantially from isothermal. And as we have 
shown previously using polytropic models of the ROSAT 
PSPC data, such small deviations from isothermality have 
a negligible impact on the derived halo ellipticities (Buote 
& Canizares 1994). 

A model X-ray image is created by integrating j x 
along the line of sight. (We ignore any spatial variations 
in the plasma emissivity convolved with the spectral re- 
sponse of the ACIS-S3 detector. Since the spectral con- 
straints indicate the gas is near isothermal such variations 
in the convolved plasma emissivity are unimportant for 
our calculations.) As in §4 we will assume the galaxy is 
viewed edge-on; i.e., we do not attempt to uncover the 
true inclination in this analysis. Because the outermost 
optical isophotes are quite flattened (cllipticity ~ 0.45) it 
is unlikely that the symmetry axis is significantly inclined 
to the line of sight. For such small inclination angles the 
halo ellipticity deduced from the X-ray analysis assuming 
i = 90° is a very good representation of the true value 
(e.g., Strimpel & Binney 1979; Buote & Tsai 1995). 

For studying the ellipticity of the dark matter halo we 
focus our attention on models where the mass is strat- 
ified on concentric, similar spheroids ("spheroidal mass 
distributions"). We focus on the following density pro- 
files having asymptotic slopes spanning the interesting 
range -2 to -4: (1) p oc (a 2 + a 2 )" 1 , (2) NFW model, 
p oc a~ 1 (a s + a)~ 2 (Navarro et al. 1997), and (3) Hern- 
quist model, p oc a~ 1 (a s + a)~ 3 (Hcrnquist 1990, 1992). 
Here a is the semi-major axis and a s is a scale length. For 
comparison we shall also summarize results for mass dis- 
tributions corresponding to potentials that are themselves 
stratified on concentric, similar spheroids ("spheroidal po- 
tentials"). Such potentials are generated by mass distri- 
butions whose ellipticity varies as a function of a; e.g., 
the logarithmic potential, $ oc ln(a 2 + a 2 ), discussed by 
Binney & Tremaine (1987). 

Given a particular mass model we compute (1) its po- 
tential, (2) p gl and then (3) the model X-ray image as 
described above. Next we compute the ellipticities and ra- 
dial profile of the model image just as done for the actual 
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Chandra data. We obtain the optimum model parameters 
by performing a \ 2 fit of the model cllipticities and radial 
profile to the data. The free parameters for the models are 
e, a s , r, and the normalization of the radial profile. Here e 
is the ellipticity of the mass for the spheroidal mass distri- 
butions or the ellipticity of the potential for the spheroidal 
potentials. 

We calculate errors on the model parameters using a 
Monte Carlo procedure similar to that done for the Geo- 
metric Test. For a given mass model (e.g., NFW) the S x (r) 
and e x values of the best-fitting model serve as the refer- 
ence values. Then we randomize each of these reference 
values according to gaussian statistics and perform fits to 
these randomized values analogously to the real data. We 
calculate the standard deviation of the best-fitting model 
parameters obtained from 20 Monte Carlo runs. We take 
these standard deviations to be the la errors. (Note that 
since the value of e x at a given a is correlated with values 
at adjacent inner a, this procedure slightly overestimates 
the errors on e x . This primarily translates to a slight over- 
estimate of the la errors for e.) 

6.2. Results 

In Table 1 we list best-fitting parameters and x 2 values 
for the "spheroidal mass distribution" dark matter models. 
Technically, such single-component mass models represent 
the total gravitating matter (stars, gas, dark matter) in 
the system. However, as we discuss below, the dark mat- 
ter sufficiently dominates the potential so that the models 
in Table 1 effectively represent the dark matter. Only e x 
values for 35" < a < 139" are included in the fits for the 
same reasons as stated for the Geometric Test (see §4.2). 

The models give consistent best-fitting values of e m 
0.35 - 0.40 with la errors of magnitude 5%-8%. The low- 
est ellipticity estimate is provided by the oblate p ~ a~ 2 
model for which the 3a error range is 0.27 < e < 0.45. 
The oblate NFW model provides the highest ellipticity es- 
timate with a 3a error range 0.35 < e < 0.47. Overall, 
though, there is very little difference in the fitted param- 
eters for oblate versus prolate models. 

Previously, in Buote & Canizares (1997) we presented 
our final analysis of the ROSAT PSPC data of NGC 720 
and inferred 0.44 < e < 0.68 (90% conf.) for the p ~ a~ 2 
model. These cllipticities are systematically higher be- 
cause the point sources discovered (and removed) by Chan- 
dra were unresolved and therefore contaminated the PSPC 
analysis. However, in that same paper we also estimated 
the effect of such unresolved discrete sources on e by as- 
suming the discrete sources would have the same spatial 
distribution as the optical light. We normalized the X- 
ray emission of such a model using the hard (T > 5 keV) 
brcmsstrahlung component discovered by ASCA. For our 
maximal discrete model we obtained 0.29 < e < 0.60 (90% 
conf.) which is fully consistent with the Chandra values. 

Recall that the average B-band isophotal ellipticity is 
(fopt) = 0.31 which is near the 3a lower limit for the 
p ~ a~ 2 models. The maximum optical ellipticity of 
« 0.45 (see Figure 3) is near the 3a upper limit for all 
the models in Table 1. Hence, the ellipticity of the dark 
matter inferred from the Chandra image is not very dif- 
ferent from that of the optical light. 

The p <~ a~ 2 model provides the best x 2 fits, and the 



quality of the fit is progressively worse for models with 
steeper density profiles; i.e., the Hernquist model fits even 
worse than the NFW model. The NFW and Hernquist 
models produce lower quality fits in part because they 
predict radially decreasing e x profiles in contrast to the 
essentially constant e x profile generated by the p <~ a~ 2 
model which agrees better with the Chandra data (Figure 
4). More importantly, the NFW and Hernquist models 
have difficulty fitting the inner bins of the radial profile. 
(The Hernquist model deviations from the data are similar 
to, but larger than, those of the NFW model shown in Fig- 
ure 4.) If the 3 inner bins are excluded the NFW model fits 
as well, yields the same e values, and predicts essentially 
the same constant e x profile, as the p ~ a~ 2 model. This 
behavior is easily understood because the best-fitting scale 
length, a s = 42", is also much larger than before which in 
effect flattens out the NFW density profile so that approx- 
imately p ~ a~ 2 over most of the image. 

Let us now consider the effect of the mass associated 
with the visible stars on the inferred ellipticity of the dark 
matter. At the very least we know that the values of e in 
Table 1 represent lower limits for the dark matter. This is 
because mass distributed like the stars predicts e x values 
smaller than observed (Figure 3). So if the stars make a 
sizable contribution to the shape of the gravitational po- 
tential, then the dark matter will have to be flatter than 
inferred in Table 1 to compensate for the rounder stellar 
potential. (Note the mass of hot gas is much less than the 
gravitating mass over the region being considered - see 
Buote & Canizares 1994). 

We also know from the Geometric Test, and from the 
results for the Hernquist models in Table 1, that a stellar 
mass model is a very poor fit to the Chandra data. Recall 
from §4.2 that the stellar mass is represented by a Hern- 
quist model with a s — 34.5" and e = 0.31. In Table 1 we 
see that even the best Hernquist models provide the worst 
fits, and since the stellar values of e and (especially) a s 
are many standard deviations away from the best-fitting 
values in the table, it is clear that the stellar model is a 
far worse fit to the data than the best p ~ a~ 2 model. It 
follows that the stars cannot make a substantial contribu- 
tion to the potential defined by the models in Table 1 as we 
found previously when modeling the ROSAT PSPC data 
(Buote & Canizares 1994). Thus, the models in Table 1 
should be associated with the dark matter. 

It is interesting to consider this conclusion in the con- 
text of the mass-to-light ratio, M/Lb- For a distance 
of 25/ifo 1 Mpc the B-band luminosity of NGC 720 is 
L B = 3.3hjQ x 10 10 L Q (Donnelly et al. 1990). Assum- 
ing a gas temperature of T = 0.6 keV (§5), the best-fitting 
oblate p ~ a~ 2 dark halo model gives M/Lb values in 
units of Ii-jqMqI Lq of 6.0 at r = R e , 12.3 at r = 2R e , and 
18.7 at r — 3R e ; here the masses are spherically averaged 
and R e — 52" = 6.2 kpc is the optical effective radius 
used to define the stellar model in §4.2. For comparison, 
the best-fitting oblate NFW model gives M/Lb values in 
solar units of 6.8 at r — R e , 12.7 at r = 2R e , and 17.0 at 
r = 3R e . 

For r — R e these values are comparable to the value 
of M/Lb ~ 7M /L determined from stellar dynamical 
studies within lR e (e.g., van der Marel 1991). The agree- 
ment of the X-ray and optical values lends further support 
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Fig. 4. — (Left panel) Radial profile denoted by diamonds (blue). Also shown are the best-fitting radial profiles generated by (oblate) DM 
halos corresponding to a p ~ a -2 profile (filled circles - red) and an NFW profile (open squares — green). (Right panel) X-ray ellipticities 
predicted by these DM models. The error bars (blue) are the values of e x measured from the source-free Chandra data, and the dashed line 
is the prediction if mass follows the stars as in Figure 3. We express a in kpc on the top axis. 



Table 1 

Ellipticity of the Dark Matter Halo 



Model 


Oblate 

a s 

e (arcsec) T x 2 


Prolate 

a s 

e (arcsec) T \ 2 


p ~ a z 

NFW 

Hernquist 


0.37 ±0.03 1.9 ±0.2 7.95 ± 0.05 32.6 
0.41 ± 0.02 26.6 ± 0.7 5.95 ± 0.03 69.8 
0.38 ± 0.03 58.0 ± 1.2 5.46 ± 0.03 99.3 


0.36 ± 0.02 2.1 ± 0.2 7.82 ± 0.05 32.4 
0.39 ± 0.02 29.3 ±0.6 5.95 ± 0.03 72.5 
0.36 ± 0.02 63.4 ± 1.3 5.49 ± 0.03 102.2 







Note. — Best-fitting values and la errors obtained for selected dark matter halo models. The mass distributions are perfect spheroids with 
ellipticity denoted by e; see §6.1 for definitions of a s and T. There are 26 degrees of freedom associated with each fit. 



to the hydrostatic equilibrium assumption. Since the X- 
ray ellipticities and radial profile (as mentioned above) in- 
dicate that the dark matter dominates the mass over most 
of the region we have studied, the agreement between the 
X-ray and stellar dynamical masses of NGC 720 also im- 
plies that the stellar dynamical mass estimate of NGC 720 
must mostly reflect the dark matter and not simply the 
stellar mass. 

Finally, the "spheroidal potential" models give results 
fully consistent with the spheroidal mass distributions. For 
example, for the model, $ oc ln(a^ + a 2 ), we obtain (for 
the oblate case) a s = 5.2" ± 0.3", T = 5.33 ± 0.05, and an 
ellipticity of the potential of 0.15 ± 0.01. Therefore, the 
best-fitting ellipticity of the dark matter is 0.20 for a <C a s 
and 0.43 for a ^> a s ; i.e., e « 0.35 — 0.40 over most of the 
image. 

7. CONCLUSIONS 

We have analyzed a new Chandra ACIS-S observation 
of the elliptical galaxy NGC 720 to verify the existence of 
a dark matter halo and to measure the flattening of the 
dark halo. The ACIS image reveals over 60 point sources 
embedded in the diffuse emission from the hot gas. Most 
of these sources were not detected with previous X-ray 
observations, and analysis of their properties is presented 
in a separate paper (Jeltema et al. 2002). We found it 



difficult to accurately replace point sources located near 
the center of the galaxy with local diffuse emission which 
hindered reliable computation of e x and PA for a < 30". 
This was not a problem for larger a where the standard 
CIAO software tool DMFILTH was able to accurately replace 
point sources with local diffuse emission. 

For 30" < a < 150" the X-ray ellipticity is consistent 
with a constant value, e x « 0.15, which is systematically 
less than the values 0.2-0.3 obtained from previous ROSAT 
PSPC and HRI observations. This discrepancy can be at- 
tributed directly to the unresolved point sources in the 
ROSAT observations. The magnitude of the PA twist 
discovered by ROSAT is confirmed by the Chandra data: 
PA w 110° for 70" < a < 150" and PA « 140° (similar 
to optical value) for a < 30". However, the twist between 
a = 30" — 60" is gradual in the source-free Chandra data 
in contrast to the sharp jump at a w 60" observed by 
the ROSAT HRI. The behavior seen by the HRI is also 
attributed directly to unresolved point sources. 

7.1. Evidence for a Triaxial Dark Matter Halo 

The X-ray PA twist for a < 150" provides immediate 
evidence for dark matter because the optical isophotes dis- 
play no evidence for such a twist. (This is a consequence 
of the "X-ray Shape Theorem" which assumes that the 
hot gas is single-phase - but not necessarily isothermal - 
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and in hydrostatic equilibrium - see §4.) The revised PA 
twist also provides strong evidence that the dark matter 
is triaxial. Previously, Romanowsky & Kochanek (1998) 
constructed triaxial models of the ROSAT PSPC and HPJ 
data of NGC 720. Their best model produced a PA twist of 
- 20° such that PA w 140° at the center and PA w 120° 
for a > 60". In Romanowsky & Kochanek's model the 
PA falls immediately, though gradually, from a = un- 
til it levels off for a > 60". This gradual twist disagreed 
strongly with the sharp jump displayed by HRI twist near 
a « 60" which the Chandra observation has shown to be 
the result of unresolved point sources in the HRI data. 

Visual comparison of the Chandra PA twist (Figure 3) to 
Romanowsky & Kochanek's best triaxial model (see their 
figure 7) reveals that they agree quite well for a < 150". 
In fact, we expect that even better agreement would be 
achieved with the Chandra data because Romanowsky & 
Kochanek had to try to match the larger ROSAT elliptici- 
ties which inhibit the ability to produce a large PA twist - 
a point we originally made in Buotc & Canizares (1996a) 
for NGC 720. Thus, the more gradual X-ray PA twist 
revealed by Chandra provides strong evidence for a triax- 
ial dark matter halo and lends important support to the 
assumption of hydrostatic equilibrium for a < 150". 

7.2. Ellipticity of the Dark Matter 

To constrain the ellipticity of the dark matter halo we 
considered both oblate and prolate spheroidal mass mod- 
els to bracket the full range of (projected) ellipticities of 
a triaxial ellipsoid. Spheroidal models have the important 
advantage over triaxial models of being less expensive to 
compute. 

First, to complement the evidence for dark matter pro- 
vided by the PA twist, we performed a "Geometric Test 
for Dark Matter" to examine whether mass distributed 
like the stars (M oc L+) can match the observed X-ray 
ellipticities (irrespective of the PA twist). Focusing on e x 
measured for semi-major axes 35" < a < 139" we find that 
the M cx L+ hypothesis is inconsistent with the Chandra 
ellipticities at the 96% confidence level for oblate models 
and at the 98% confidence level for prolate models. Thus, 
both the PA twist and the ellipticities of the Chandra im- 
age imply the existence of dark matter independent of the 
temperature profile of the gas. This evidence for dark mat- 
ter cannot be explained by general modified gravity theories 
such as MOND. 

Second, for an assumed spheroidal model of the dark 
matter we generated a model X-ray image to compare the 
model ellipticities and radial profile to the Chandra data. 
The dark matter density model, p oc (a 2 + <z ) , provides 
the best fit to the data and gives ellipticities and la errors 
of e = 0.37 ± 0.03 (oblate) and e = 0.36 ± 0.02 (prolate). 
NFW and Hcrnquist models give similar ellipticities for 
the dark matter. These ellipticities are less than the val- 
ues of 0.44-0.68 (90% conf.) obtained from our previous 
ROSAT observations which did not account for unresolved 
point sources. In our previous ROSAT and ASCA analysis 
of NGC 720 which did attempt to account for unresolved 
point sources assuming such sources follow the stellar light 
(Buote & Canizares 1997), we obtained 0.29 < e < 0.60 
(90% conf.) which is fully consistent with the (source- 
free) Chandra values. It follows that unresolved sources 



remaining in the Chandra image should have a negligible 
impact on the dark matter ellipticities determined in this 
paper. (This conclusion can be checked directly with a 
deeper Chandra observation so that point sources can be 
detected to a lower flux limit.) 

The ellipticity of the dark matter inferred from the 
Chandra data is not very different from the optical light, 
(e pt} = 0.31 and e™ a t x w 0.45. But this does not con- 
tradict the results from the Geometric Test. The primary 
reason why the MccL 4 model cannot produce the X-ray 
ellipticities is because L± is too centrally concentrated. 
Consequently, the spherically symmetric monopole term 
dominates the stellar potential and predicts values of e x 
that are smaller than observed, particularly for a > 60". 

The NFW and Hcrnquist mass models do not fit nearly 
as well as the p <~ a~ 2 model, this discrepancy primar- 
ily arising from the radial profile in the inner regions. 
If the three inner bins are excluded from the fits (i.e., 
a < 16") then the NFW model fits as well as the p <~ a~ 2 
model. But the resulting larger scale-length of the NFW 
model implies that the NFW density behaves approxi- 
mately as p ~ a~ 2 over the a range fitted. So unless 
the mass actually consists of two components, i.e., "some- 
thing" plus NFW, the p ~ a~ 2 model is clearly favored. 
This "something" would have to be more centrally con- 
centrated than the optical light (which follows an i? 1 / 4 
law with R e = 52"). 

7.3. Strange Behavior at Large Radius 

At large radius (150" < a < 180") e x increases steadily 
to a value of « 0.20 while the position angle decreases 
to a value of « 80°. These changes are highly statisti- 
cally significant, although (1) the fact that this region is 
near the edge of the CCD, and (2) the PAs measured by 
ROSAT and Chandra in this region do not agree, suggest 
that an independent observation with the wide-field CCDs 
on XMM should be used to verify the results in this region. 

If the Chandra values for e x and PA are correct in this 
region, then a physical interpretation is necessary. We 
have argued that there are several good reasons to expect 
that the hot gas is in hydrostatic equilibrium (see below). 
However, it is possible that the hydrostatic assumption 
is valid only for a < 150" where the PA twist is consis- 
tent with the triaxial dark matter halo of Romanowsky & 
Kochanek (1998); i.e., the divergence in PA for larger a 
could signal the departure from hydrostatic equilibrium. 

If the hydrostatic condition does hold for a > 150" then 
there must be another, larger scale, component in the dark 
matter. The possibility of large-scale dark matter in NGC 
720 has been suggested from analysis of its dwarf satel- 
lites by Dressier et al. (1986). Why such a separate dark 
component decides to assert itself at a « 150" is unclear. 

7.4. Hydrostatic Equilibrium vs. Merger 

Since the evidence for dark matter in NGC 720 pro- 
vided by the Chandra data rests primarily on the assump- 
tion that the hot gas is in hydrostatic equilibrium, here 
we summarize the arguments in favor of equilibrium: (1) 
external influences on the diffuse gas should be negliblc be- 
cause NGC 720 is very isolated (§4.1), (2) other than the 
PA twisting, there are no notable asymmetries in the X- 
ray isophotes - in particular no substantial centroid shifts 



13 



(§4.1), (3) an N-body / hydrodynamical simulation sug- 
gests for a system without obvious subclustering the X-ray 
shape analysis should be valid (§4.1), (4) evidence for dark 
matter is provided separately by e x and PA (§7.1 and §7.2), 
(5) Romanowsky & Kochanek triaxial model explains PA 
twist within a ~ 150" (§7.1), and (6) the gravitating mass 
within li? e determined from X-rays agrees with the stellar 
dynamical value (§6.2). 

A final argument for equilibrium, at least in the central 
regions, is provided by the ages of the stars. Terlevich & 
Forbes (2002) estimate a relatively young age of 3-4 Gyr 
for the central stellar population in NGC 720, consistent 
with a fairly recent merger with a gas-rich neighbor. How- 
ever, the sound crossing time of the hot gas in NGC 720 
is only w 0.1 Gyr indicating that the gas has had ample 
time to settle down since the starburst. (Also, in the opti- 
cal NGC 720 is classified as a "core galaxy" (Lauer et al. 
1995) and has rounder isophotes near the center suggesting 
that it has not been dynamically disturbed in a long time 
according to the simulations of Ryden & Tinker (2002).) 
We conclude that for a < 150" (18.2h^ kpc) that hydro- 
static equilibrium is very likely to hold sufficiently well for 
analysis of its gravitating matter distribution. 

The only evidence that may suggest a recent merger is 
associated with properties at large radius. First, the PA 
twist at large radius (a > 150") cannot be explained by a 
single triaxial ellipsoid in hydrostatic equilibrium (§7.3).. 
Second, the "arc" of point sources at a ~ 120" in NGC 
720 discussed by Jeltema et al. (2002) is similar to shells 
seen in some elliptical galaxies which are thought to be 
merger remnants. 

7.5. Implications for Galaxy Formation and Cosmology 

The existence of dark matter in isolated elliptical galax- 
ies has been assumed by the standard Cold Dark Matter 
(CDM) paradigm for many years, but primarily because of 
the ambiguities associated with the analysis of stellar ve- 
locities, definitive evidence for dark matter has been lack- 
ing. The Chandra observation of NGC 720 provides vital 
evidence that dark matter does indeed exist in an isolated 



elliptical galaxy. This geometric evidence for dark matter 
is apparently the first dynamical tracer in any stellar sys- 
tem that has successfully distinguished dark matter from 
MOND (cf. Scllwood & Kosowsky 2002a). 

The evidence for flattened dark matter provided by the 
cllipticities of the Chandra X-ray image of NGC 720 is in- 
consistent with the Self-Interacting Dark Matter (SIDM) 
model of Spergel & Steinhardt (2000). The SIDM model 
predicts that the dark matter halo should be spherical 
within a radius of ~ 10 kpc in NGC 720 using the "colli- 
sional radius" defined by Miralda-Escude (2002). On the 
other hand, the moderate flattening of the dark matter 
halo in NGC 720 is inconsistent with the highly flattened 
distributions predicted by the cold molecular dark matter 
model (Pfenniger et al. 1994). 

The inferred cllipticity e « 0.35 — 0.4 of the dark matter 
halo of NGC 720 is consistent with simulations of galaxy- 
scale halos in a ACDM model, but may be larger than 
indicated by a AWDM (i.e., "warm" dark matter) model 
(Bullock 2002). These ellipticities also are similar to the 
small number of measurements that have been made in 
other types of galaxies using different methods (e.g., Mer- 
rifield 2002). Finally, the evidence for triaxiality provided 
by the X-ray PA twist is another important verification of 
CDM models (e.g., Bullock 2002). 

Although the Chandra observation of NGC 720 has pro- 
vided important evidence for flattened, triaxial dark mat- 
ter in an elliptical galaxy it is essential to apply this anal- 
ysis of X-ray shapes to other galaxies (and groups and 
clusters) to determine whether the conclusions drawn from 
NGC 720 are representative of other systems. Given the 
importance of detecting and removing point sources for 
the analysis of dark matter in NGC 720, it is clear that 
Chandra and its « 1" resolution will be essential for such 
studies of other galaxies. 

We thank M. Krauss for helpful discussions regarding 
the features in the ACIS-S exposure map, and J. Sell- 
wood and the anonymous referee for comments on the 
manuscript. TEJ was supported by an NSF fellowship. 
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